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Abstract 

We study at LEP II and Next Linear Collider energies Higgs production via the 
bremsstrahlung channel e"*'e~ ZH, with Z — > /i^/i^ and H bb, and the cor- 
responding irreducible background, in presence of hard photon radiation, both from 
the initial and the final state. We carry out an analysis that includes the compu- 
tation of all the relevant contributions to the complete tree^evel matrix element 
for e'^e~ — > fi^fi~bbj and makes use of the one at leading order e^e~ — > ^jl^ ^~bb 
interfaced with electron structure functions. We concentrate on the case of mass 
degeneracy Mh ~ Mz, for which next -to-leading electromagnetic contributions can 
modify the content of 66-pairs coming from H and Z decays. A brief discussion 
concerning the case Mh ^ Mz is also given. 



^E-mails: Moretti@to.infn.it; Moretti@hep.phy.cam.ac.uk. 



1. Introduction 

If the Standard Model (SAi) Higgs boson has a mass around 90 GeV, it is degenerate with 
the neutral electroweak (EW) boson Z {Mz ~ 91.175 GeV). Although the model itself 
does not show any preference for this value, such a condition could well occur in nature. At 
this mass the Higgs particle predominantly decays into 66-pairs, and this is also a possible 
signature of the Z-vector. Therefore, whenever the two bosons are contemporaneously 
produced, their peaks overlap in the invariant mass spectrum of the decay products and 
cannot be separated, not even for optimistic detector performances^. However, a closer 
inspection of the candidate 66-sample (in a window around Mh ~ Mz) should reveal a 
6-quark content larger than the one arising from Z-decays only. This would then allow 
one to recognise the presence of the Higgs boson. To this end, the fact that the Higgs 
Branching Ratio (BR) into 6-quarks is ~ 90% (for Mh ~ 90 GeV), whereas Z ^ bb 
constitutes only 22% of the possible hadronic Z-decays, is very helpful^. 

A Higgs boson with Mh ~ Mz was beyond the discovery potential of LEP I. Unless the 
Tevatron is able to detected it in the short term, the first machine where such a particle 
can be searched for is LEP II, provided that its Centre-of-Mass (CM) energy reaches 
~ 200 GeV. In fact, at this energy, the machine will be beyond the threshold of an on- 
shell ZiJ-pair, which is the dominant Higgs production mechanism at LEP II. Recently, 
it has been assessed ^ that the cleanest Higgs signature at LEP II is the one proceeding 
through the decays Z —>■ ^~ and H ^ bb (the 'golden channel'), although also other 
final states can be considered {H t^t^ and Z e"'"e~, i/z/, jets). In the more distant 
future, high precision Higgs physics can be performed at the Next Linear Collider (NLC, 
with a/s = 350 — 500 GeV) [0. An e^e~ linear machine has a clear advantage with respect 
to hadron colliders. Because of the reduced QCD background at the former, one can more 
easily exploit the H ^ bb decay channel, whereas at the latter extremely high 6-tagging 
efficiency and not-6-jet rejection are required in order to make Higgs detection feasible^. 

At LEP II the case Mh ~ 90 GeV will require the highest CM energy option, which 

^In fact, the Higgs width Th, for Mh = 90 GeV, is w 2.4 MeV, whereas experimental resolutions are 
generally comparable or larger than the Z width, Tz = 2.5 GeV. 

•^An experimental method of discriminating between jets that originate from &-quarks and those pro- 
duced by light quarks and gluons, namely fo-tagging, has been successfully exploited in the past few years, 
especially at LEP I and the Tevatron. The method requires one to tag hadronic events with secondary 
vertices. 

''At the moment, it is not clear whether this will be possible or not, particularly in case of Z-H 
mass degeneracy. For discussions about the Higgs discovery potential of future hadronic machines, see 
Refs. |,|. 
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corresponds to a second stage of the project of the machine, therefore quite far in the 
future. Furthermore, by assuming an yearly luminosity of 500 pb^^ for the CERN collider 
and of 20 fb^^ for the linear machine, one ends up with a number of events that is larger 
in the second case. Also, whereas the gauge and Higgs bosons produced at LEP II would 
be practically at rest, such that their decay products are generally soft and uniformly 
distributed over the accessible phase space, at NLC energies these have a strong boost, 
making their tagging relatively easier. For all these reasons we believe that the NLC will 
still play a decisive role in detecting and, particularly, in studying the properties of a Higgs 
boson with a mass near Mz- 

We earlier mentioned that the Zif-channel is the largest Higgs production mechanism 
at LEP II. If Mh ^ Mzi at ^/s ~ 350 GeV it is still the dominant process, whereas at y/s ~ 
500 lyiy-fusion has bigger rates (ZZ-fusion is smaller than e^e" —>■ ZH) p. However, 
the Zif-mechanism has a few advantages with respect to the WW-one: its final state 
is fully constrained and it allows Higgs spin, parity and branching ratio determinations 
(see 1^, and references therein). Thus, the ZH-process is of a crucial importance even 
in energy regimes where it is not the largest Higgs production mechanism. This is the 
reason why we concentrate here on the bremsstrahlung process only. In addition, we 
restrict our attention to the case H ^ bb (as the corresponding BR is much larger than 
that one into taus). Concerning the Z-boson, we select muon decays only, thus avoiding 
complications due to hadronic final states (when Z goes into jets or taus), to additional 
t-channel backgrounds (for Z —>■ e^e~) and to the loss of kinematical constraints (final 
state neutrino). 

The story of theoretical studies of the process e~^e^ ZH dates back to Refs. [|, |, 
where the Born cross section was first computed. In Ref. the leading order (LO) cross 
section for e~^e~ — >■ ffH was calculated analytically. Next-to-leading order (NLO) EW 
corrections were computed later 1T|, |T|, |16|, As for Mh ~ 90 GeV the 

Higgs width is a few MeV only, the Higgs boson can be safely kept on-shell in the 
computations. However, for phenomenological studies, the full final state, including the 
Higgs decay products, has to be known. Therefore, signal and background processes of 
the type, e.g., e~^e~ fi^fi'bb, must be considered. Recently, tree-level semi-analytical 
cross sections for e+e~ — ^ fi'^fi~bb events were given [|1^, Alternatively, one can 

exploit a Monte Carlo ^ or a 'deterministic' approach p2 . 



It is the purpose of this paper to quantify the influence on the cross section of hard 



^For a review, see (and references therein). 
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photon emission, which can take place in the process e'^e^ fi^fi^bb. In fact, contrary 
to LEP I, where the Z-resonance imposes a natural cut-off on events with hard photons 
produced by the initial state, at LEP II and NLC such a suppression does not act any longer 
and, in addition, as the beam energy is larger, the probability that the incoming electrons 
and positrons can radiate increases. It follows that a sample of pure e+e^ fi^^^bb 
events, without 7-radiation, does not exist in practise and one has to deal inevitably with 
electromagnetic (e.m.) emission. We also notice how all particles in the final state are 
charged and therefore can also radiate hard photons. 

The inclusion of higher order e.m. effects is especially important in the case ^ Mz, 
since radiative corrections modify the relative content of H and Z-decays in the candidate 
66-sample. This is true whether one adopts missing mass techniques by computing the 
invariant mass recoiling against the tagged /i"''/i~-pair, or one attempts the reconstruction 
of the 66-signal directly from the 6-jets p3|. In the first case, photons enter in the Mrecoii 
spectrum and also spoil the reconstruction of the Z-boson via the tagged muons, since 
these latter can emit e.m. radiation after the Z-decay. In the second case, as 7-radiation 
can take place also off 6-quark lines in the Higgs decay process, such effects are equally 
present (although with a reduced probability since \Qb\ = 1/3 < \Qe,fi\ = 1 and because 
only half of the diagrams occurring in case of radiation from leptons are involved here). 
To avoid this, one could try to include photons by looking at the Mf,i^ invariant mass. 
But such a procedure would also be distorted by unwanted 7-radiation produced by e+e" 
and pairs. Moreover, the knowledge of the absolute size of the corrections for real 

photon emission could help if one discriminates between a SAi Higgs boson and a possible 
neutral scalar from Supersymmetric models, and/or in testing the Yukawa coupling of the 
Higgs particle. 

In general, the hard e.m. radiation affects the kinematics of the ^^fi^bb events in vari- 
ous ways. First, the presence of a hard photon modifies the two-body decay kinematics of 
the process e~^e^ ZH fi~^fi~bb ||2^, as photon emission cannot be unambiguously 
assigned to one or the other of its three possible sources: e~^e~ , jX^n' and bb. Second, the 
radiated photon will affect a cut around Mz on the /i"'"/i~-pair, which allows one to get rid 
of a large part of the backgrounds in the non-radiative case. In fact, many of the config- 
urations coming from the ZH signal do have muons whose invariant mass does not peak 
at the Z-pole. In this respect, as the radiation from the e-lines is concentrated to a large 
extent along the beam direction and that from 6-quark lines is suppressed by the charge 
Qbi as a consequence, the /i-lines are the main source of detected photons. Finally, the 
number of background diagrams leading to fi'^fi~blry final states in e~^e~ initiated processes 
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is very large (168 Feynman graphs at tree-level, neglecting 62 diagrams which include a 
suppressed if/i/i Yukawa coupling), compared to the number of those involving the Higgs 
resonance (6 graphs at tree-level)^. 

In carrying out our analysis, we will make some simplifications, which should not modify 
the conclusions we will get in the end. We will perform only a tree-level calculation of 
the processes e"^e^ —>■ fi^fi^bb and e~^e~ —>■ fi^fi~'blry (with the photon emitted both from 
the initial and the final state). We will however 'dress' the former with the Initial State 
Radiation (ISR) |2^, implementing this latter by means of a convolution with electron 
structure functions (in particular, we use the expressions given in Ref. |^)[|. This is 
done in order to be able to sum up the rates corresponding to the two above processes, 
through the order 0{al^). In fact, the ISR formulae, as given in Ref. |]28|, include the exact 
photon corrections to the e'^e~ annihilation subprocess up to the second order in the e .m. 
coupling constant (in particular, they also embody hard photon emission). For consistency, 
we implemented only 0{aem) ISR terms (the inclusion of the C(ag^) pieces does not 
change significantly the size of such corrections). By computing the total cross section 
of the leading (2 — 4) reaction supplemented by the ISR and that one of the radiative 
(2 — > 5) process, we are then able to separate at the order the contribution due to 
hard photon emission from the one due to collinear and soft e . m . radiation (see below 
the cut in transverse momentum) and to virtual corrections. In particular, we combine 
the two approaches by subtracting the radiative (2 — > 5) Matrix Element (ME) from the 
leading log 0{aem) part of the ISR, to avoid double counting (see also Ref. |^ 0)- In the 
following, the label 'LO' will identify rates obtained from the 2 — > 4 process dressed with 
ISR and with the hard emission subtracted, whereas 'NLO' will refer to rates obtained 
from the 2 — 5 reaction involving hard e . m . radiation only. This subtraction is performed 
between, on the one hand, the graphs in Fig. 1 and, on the other hand, those in Figs. 2-4 
in which photons are connected to electron and positron lines: i.e., 1 & 2 in Fig. 2, 1-2, 
7-9 & 14 in Fig. 3, 1-2 & 8-9 in Fig. 4. Since in neutral current processes ISR and Final 
State Radiation (FSR) are separately gauge invariant, this procedure is legitimate here. 

The above procedure is adopted for the case of the dominant, universal, factorisable and 
process independent contribution to the complete set of ISR corrections to four-fermion 



^At leading order in e^e — > /i^/i bb events one has 24 and 1 diagrams, respectively, again excluding 
9 diagrams with H^^ couplings. 

^We neglect to consider Linac energy spread and beamsstrahlung, as they are negligible compared to 
the ISR, at least for the 'narrow' D-D and TESLA collider designs These effects would however be 
straightforward to insert. 
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production, containing all mass singularities L = ln(s/mg) and expressed via the ISR 
'radiator', as known from s-cliannel e^e~ annihilation^. In addition to these, there are 
also non-universal, non-factorisable and process dependent ISR corrections, which arise 
in association with t- and w-channel electron exchanges (neutrino, in case of charged 
currents). The full set of formulae needed to incorporate the complete ISR in CCS, NC2 
and NC8 e+e~ — > 4 fermion processes has been recently given, in Ref. |3^. Whereas 
Higgs signals (graph d in Fig. 1) and non resonant background contributions (graphs b and 
c in Fig. 1) proceed via s-channel annihilations, the remaining background contributions 
(graphs a in Fig. 1) do through t- and w-electron channels, involving ZZ-, Z7- and 
77-production (i.e., NC8-type diagrams). Thus, we could well apply non-universal ISR 
corrections to these latter contributions. However, it has been shown in Ref. |^2| that 
they are generally smaller by an order of magnitude with respect to the universal ones. 
Furthermore, one has to consider that we are concerned with Higgs searches for Mh ~ Mz, 
such that in the end (after the Higgs selection procedure) the dominant background is ZZ- 
production (i.e., NC2-type graphs), for which the non-universal ISR corrections are even 
more suppressed. For these reason then, we neglect here such effects. 

Compared to the universal ISR, QED corrections related to the final state (FSR) are 
much smaller (and comparable to the non-universal ISR corrections I^Q], |3^, |3^). 



Typical suppressions of the order 0{r/M) are expected for the interferences between 
ISR and FSR in resonant boson pair production as well as in inter-bosonic final state 



interferences We do not expect these effects to be relevant for Higgs production 

either. Thus, because of their small size, we do not include them in our analysis. 

We also ignore genuine weak corrections, as these have been proved to be well under 
control (at least for the on-shell process e^e~ —* ZH) [p6| . 

We further stress that, on the one hand, a calculation including even part of the above 
effects would be extremely CPU-time consuming, because of the large number of Feynman 
diagrams involved in the processes studied, which have up to seven external particles and 
multiple Breit-Wigner resonances in different regions of phase space. On the other hand, 
a full 0{alj^) calculation is well beyond the intentions of this paper. 

A further simplification we have adopted is to avoid computing some of the interference 

terms of the two processes e~^e~ — ^ fi~^fi~bb{''^), which either vanish identically or are 

extremely small compared to the squared terms of the MEs. Finally, we do not discuss 

the QCD background, e~^e~ Z + n jets (with n > 2), as this is largely suppressed if all 

^Naively, ISR universal effects tend to reduce the effective beam energy, thus modifying both the 
normaUsation and the shape of the differential distributions which are of interest in Higgs searches. 
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jets are well separated (as we assume throughout this paper) and they have to reproduce 
the Higgs mass, Mjets ~ Mh- 

Finite width effects for the Z and H are included|. The masses of both 6-quarks and 
muons are non-zero, apart from the case of the H^^ coupling, in which has been 
set equal to zero, thus eliminating the corresponding diagrams from the MEs. No effort 
has been made to simulate experimental effects, in particular in detecting 6-quarks, apart 
from assigning plausible values for the 6-tagging efficiency and the misidentification 
probability tc of c as a h. We cut on the transverse momentum (with respect to the beam 
direction) of each particle in the final state, pfJ^ > 1 GeV, which, on the one hand, is 
reasonably compatible with detector requirements and, on the other hand, protects our 
radiative MEs from both divergences and numerical instabilities. Hence, this value of 
transverse momentum also defines the phase space region of soft and collinear photons 
(i.e., < 1) emitted from the incoming e"'"e~-lines where we make use of the ISR (as 
described above). 

In the next Section we give an account of the numerical calculations; in Section 3 we 
present and discuss our results and Section 4 gives our conclusions. 

2. Calculation 

For Mh ~ both at LEP II and at the NLC the cleanest Higgs signature in the 
intermediate mass range will be via the process 

e^e~ fi^fi~bb. (1) 

The corresponding Feynman diagrams for background (graphs a, b and c) and signal (graph 
d) contributions are depicted in Fig. 10. Among the background graphs, one distinguishes 
between 'conversion diagrams' (Fig. la) and 'annihilation diagrams' (Fig. lb,c). They are 
also called crab and deer diagrams, respectively, further dividing these latter into /i-deer 
(Fig. lb) and 6-deer (Fig. Ic) diagrams |TB|, By inserting the appropriate type of 

internal propagator (7 or Z) and by performing all the possible permutations of the boson 
lines on the fermion ones one can get out of the background diagrams of Figs, la-c a total 
of 24 diagrams, 8 per each kind of graphs (a, b and c). The signal diagram in Fig. Id is 
unique. 

^We are not concerned about possible gauge invariance violations due to bremsstrahlung off unstable 
particles, as there is no cancellation here which could amplify such effects. 

^°Here and in the following we adopt the labelling: e+ (1), e~ (2), (3), /i^ (4), b (5), b (6) and 7 

(7). 
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To obtain the Feynman graphs representative of the radiative process 



e^e — >• bb'j, (2) 

one has to add a real photon in all possible ways to the diagrams of Fig. 1. This is a 
topologically trivial matter, however, the total number of graphs to be computed is very 
large. The case of the signal is the simplest, as one obtains the 6 diagrams of Fig. 2 (here 
the wavy internal lines represent a Z). For the case of the crabs one gets 56 diagrams 
in total, as indicated by Fig. 3 (in which an internal wavy line represents both a 7 and a 
Z). This is also the number that one gets from the deers (both for fi and b, see Fig. 4). 
In total, one ends up with 174 diagrams at tree-level. As already mentioned, we did not 
consider 62 diagrams involving Hfifi Yukawa couplings. 

More than from the large number of graphs, complications in computing the cross 
section for process (H) arise from the fact that they have in general rather different resonant 
structures over the accessible phase space. For example (see Fig. 2), the signal diagrams 
present Breit-Wigner peaks of the type H ^ bb, H ^ bb'-f, Z yU/i and Z fifi'y. For 
the backgrounds, in addition to the two already mentioned Z-resonances, one also finds 
Z ^ bb and Z —>■ bb'-f peaks, together with diagrams which are not resonant at all (for 
example, the ones in Fig. 3 when all the internal wavy lines are identified with virtual 
photons)!^. The (numerical) integration must carefully take into account this diversified 
structure. 

We computed the cross sections by using the packages MadGraph and HELAS |^ 



for the Feynman diagram evaluation and VEGAS for the integration over the phase 
space. For VEGAS, we split the total Feynman amplitude squared into a sum of non- 
gauge-invariant terms (as already done in Refs. PH M), each of which has a particular 



resonant structure, and integrate them by using appropriate phase spaces which absorb 
the Breit-Wigner peaks in the integrand. In general, the change of variable 

Q^-M^ = MT tan 6, dQ' = ^-^ Jf "^^^ 

where Q, M, T stand for the virtuality, the mass and the width of the resonance, gives 

an integrand which depends smoothly on 6. For example, in Fig. 2, it is convenient to 

isolate three terms, i.e., the amplitude squared of the graphs: i) 1 & 2; ii) 3 & 4; iii) 5 & 

^^As we will be eventually interested in studying the distributions in the invariant masses Af^^ and 
Mf£^, for the time being we do not discuss here the fact that also resonances of the type Z ji^ ji^hb 
can appear in the case of the deers, see, e.g., graph 1 of Fig. 4 (in which the wavy line attached to the 
incoming fcrmions represents a Z and the one joining ^- and 6~lines is a photon). 
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6, and to compute them separately by means of different phase spaces. The interferences 
between the set of diagrams in i), ii) and iii) are the last contribution to the part of the 
total ME which includes Higgs resonant graphs. In a similar way, one proceeds with the 
background graphs in Figs. 3-4. 

Such a procedure gives a substantial reduction of the integration errors, although some 
diagrams need to be computed twice or more (e.g., diagrams 1 & 2 in Fig. 2 are necessary 
for the I Ml + term as well as for some of the interferences). This apparently increases 
the final CPU time of the run. However, the competition between these two aspects (i.e., 
high precision but multiple Feynman diagram evaluations) is such that the computing time 
needed to get a given accuracy is much less than that required to integrate the differential 
cross section without any special care, especially when one uses adaptive algorithms for 
the multi-dimensional integrations, as we have here. 

Among the different non-gauge-invariant terms one has to deal with, the interferences 
between the various sets (of signal, crabs, fi- and 6-deers) as well as the ones between 
graphs within the same set are the most complicated, as these generally mix up graphs with 
different resonant structures^. However, this often implies that they are small compared 
to the amplitudes squared of diagrams with identical composition of Breit-Wigner peaks, 
since the phase space region in which one or more graphs are large is different. Therefore, 
a useful and time-saving procedure is to neglect such terms, whenever possible. 



For example, concerning the leading order process (ffl), it has been shown in Ref. [20 



that the various Higgs-background interferences either vanish identically or are small, so 
that they can be safely neglected. In the case of the interferences between the various 



sources of background such a suppression is less visible |jT8[ . We have checked whether this 

remains true in presence of hard photon radiation. For example, by interfering the signal 

diagrams of Fig. 2 with the crabs and deers of Figs. 3 and 4 together, one obtains mixed 

terms which are negligible compared to all the squared contributions, whereas mixing the 

various sets of crab and deer diagrams of Figs. 3 and 4 yields interferences which are of 

the same order as some of the amplitudes squared. This also happens for the interferences 

within the sets of Figs. 2, 3 and 4, if treated separately: that is, when one splits the modulus 

squared of the signal process represented in Fig. 2 as previously described, and those of 

the backgrounds in Figs. 3 and 4 by isolating the Z-resonances out of the graphs with 

7-propagators. In this respect, we found it useful to recognise within the crab diagrams 

of Fig. 3 the components with the following resonant structures: i) Z ^ fi'^fi^{'y) and 

-"^^We have integrated them by using a flat phase space, which does not map any of the possible peaks 
of the interfering graphs. 
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Z — > bb{'y); ii) Z ji^ji^ij) and 7 ^{l)] iii) 7 ^ S'lid Z ^ bb{j); iv) 

7 ^ fi'^fi~{'y) and 7 — bb{'y). In the following, we will refer to them as ZZ, Z7, 7Z 
and 77 backgrounds. A Z-resonance is present also in the diagrams of Fig. 4, however, 
we will refer to these collectively as /i- and 6-deers (although in the computations their 
resonant structure has been appropriately taken into account), as these diagrams are very 
suppressed if compared to those in Fig. 3 and also because they do not have the two- 
to-two and two-to-three kinematics typical of the signal (compare to Fig. 2). As already 
stressed, we recall here that all these amplitude contributions we have been discussing must 
be summed up together in the end. However, if taken separately, they provide a useful 
way of looking inside the process and distinguishing between the different fundamental 
interactions. In summary, in the results we will present in the next section, we have 
systematically neglected interferences between the signal and the backgrounds whereas we 
have kept all the others. 

The following numerical values of the parameters have been adopted: Mz = 91.175 
GeV, Tz = 2.5 GeV, Mw = 80.23 GeV, Tw = 2.2 GeV, and for the Weinberg angle we 
have used its leptonic effective value of smlQ{9w) = 0.2320. For the fermions: = 0.105 
GeV and = 4.25 GeV. The e .m. coupling constant has been set equal to 1/128. For 
the Higgs width Th we have adopted the tree-level expression corrected for the running of 
the quark masses in the vertices Hqq (these have been evaluated at the scale /i = Mh PI ). 
Therefore, in order to be consistent, we have used a running 6-mass in the Hbb vertex 
of the production process here considered. As representative values of the CM energy of 
LEP II and NLC we have adopted 200, and 350, 500 GeV, respectively, whereas the Higgs 
mass has been fixed at 90 GeV. 



3. Results 

In the non-radiative process e^e^ —>■ fi^fi^bb Higgs signals should appear as narrow reso- 
nances in the invariant mass spectrum of the Higgs decay products over the contribution 
due to the background processes, which are generally flat, apart from the region around 
Mz, where the Z peak (due to the ZZ crab diagrams of Fig. la) clearly sticks out. In 
order to detect Higgs signals via the two-6-two-/i channel one can adopt two different 
strategies 

^■^We ignore the full reconstruction of the reaction e+e^ ZH via the decays into jets and/or lep- 
tons of both the H- and the Z-boson, by fitting their kinematics, because it would present additional 
complications, which are beyond the intentions of this study |p3| . 
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the calculation of the missing mass recoiling against the reconstructed Z — > 



by plotting the distribution in Mrecoii = \j\{Ve+ + Pe-) — + Pa*-)]^ aiid exploiting 
tagging on the yU^/i~-system only; 

• direct reconstruction of the if-peak from the 6-jets, by plotting the distribution in 

Once a hard photon is included via the process e^e^ — » ^^^^bb~f, the next-to-leading 
order contribution to the first spectrum coincides with the invariant mass of the 667- 
system, thus the rates up to the order in the recoiling mass are the sum of the 
distributions M^^ at LO and M^^^ at NLO. In the second approach, the full rates up 
to Ogj^ in the M^j spectrum are given by the sum of the 66-invariant masses at leading 
and next-to-leading order. We expect the possible alternative strategy we outlined in the 
Introduction (i.e., to look at the 'pure' M^b^ spectrum) to be less successful, because in this 
case one would be able to correctly reconstruct Higgs peaks only in the case of diagrams 
5 and 6 (of Fig. 2), which give a suppressed contribution to the total e+e~ fi~^fi~bb'-)' 
cross section. In presenting our results, we concentrate then on the two spectra M^^i and 
Mfcs^, with Mh^Mz. 

In order to suppress the backgrounds jZ and 77 (among the crabs) as well as the 
contributions from the /i- and 6-deers, we require that the invariant mass of the muon 
pair reproduces a Z-boson, by applying the cut, e.g., AMz = |M^+^- — Mz\ < 10 GeV 
P, |2^. In this way, we mainly select the ZZ contribution from the crabs, since the Z'j 
one will be largely suppressed in the end by a cut around Mh. In this way, one is able 
to compare Z ^ bb and H ^ bb decays [^. Furthermore, in order to reduce the content 
of Z-decays in the candidate 66-sample we consider events for which |cos6'^+^-| < 0.8 
[p3| (where 6*^+^- is the angle of the reconstructed Z-boson with respect to the beam), 
since e~^e~ ZZ events are strongly peaked in the forward/backward direction due to 
the t, u-channel exchange of electrons [^[5 

Tables I-III show the cross sections in a window of 20 GeV, centered around Mh = 90 
GeV (i.e., we consider events with AMh = 1^^5(7) — Mh\ < 10 GeV), for the signals 
e~^e~ ZH fi^ fi~bb{'y) {S, the square of the last diagram in Fig. 1 at LO and of those 
in Fig. 2 at NLO) and the total backgrounds e~^e~' ZH fi'^ fi~bb{'~f) {B, the square of 
the sum of the first three sets of graphs in Fig. 2 at LO and of those in Figs. 3 and 4 at NLO), 
together with the total significance S, for the integrated luminosity C = J Ldt = 0.5(20) 
fb-i per annum at LEP II(NLC), for = 200 (Table I), 350 (Table II) and 500 GeV 

^''This is also true for the cases Zj, 7Z and 77. 
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(Table III). All the mentioned constraints have been applied: that is, pj.'* > 1 GeV, 
/S.Mz < 10 GeV and | cos 6*^+^-1 < 0.8. In case of cross sections {S and B columns), the 
label NLO identifies rates from the radiative process e+e^ fi^fi^bb'-f only (with pj, > 1 
GeV). In particular, the second line represents rates obtained from the distribution in 
invariant mass of the system 667, whereas the other refers to rates obtained from the Mj,^ 
spectrum. In case of significances (S columns), the label NLO give rates obtained through 
the order (9(a;g„) (all photons) in case of the direct reconstruction method (upper line) 
and of the missing mass analysis (lower line), as obtained by adopting the subtraction 
procedure described in the Introduction. According to this treatment, LO refers here to 
rates from the process e"'"e^ — > fi^fi^bb corrected for universal ISR effects due to virtual 
and real photons (these latter with pj,<l GeV). For reference, we also give (in square 
brackets) cross sections and significances when no radiation is present, as obtained by 
integrating the spectra in Fig. 5. 

In computing the S's (see Ref. [0), we neglect the (small) probability of misidentifi- 
cation of light quarks and gluons as 6's and we assume BR{H bb) » BR{H — > cc). 
Then, the probability of picking one b out of two is [1 — (1 — e^)^], whereas the total 
significance E is 

^Ja{e+e- ^ ZH ^ ^i+fi-bb{-f)) 

with 

- ^ - - (5) 

[l-(l-6,)2]+5[l-(l-6,')2]' 

where the factor 6 has been introduced in order to remind the reader that the EW couplings 

7*,Z — i> qq{'y) are generally different, depending whether q = c or b. However, to use 

5 = 1 in formula (|^) is always a good approximation for our analysis. In fact, on the 

one hand, the background B practically coincides at tree-level with the ZZ-piece and 

BR{Z bb) ^ BR{Z —* cc). Whereas, on the other hand, the NLO rates (for which 

BR{Z — > 667) ^1/4 BR{Z — > CC7), although this does not occur in all diagrams) have 

to be added to the LO ones to produce the correct significances, such that the effect of 

the different e . m . coupling is largely washed out in the endQ. We consider the following 

five combinations of 6-tagging efficiency (ef,) and c —>■ b misidentification [ej)'- £6 = 1 and 

ec' = (perfect tagging), and e^' = 0.2eb [H, for et, = 0.2, 0.4, 0.6, 0.8. 

^^The value of 6 would be different from 1 when using the NLO rates on their own, but E's computed 
in this way have no meaning for the present analysis. 
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Before proceeding further, a few points concerning the background due to top events are 
worth mentioning here. In fact, tt-production followed by the decay ti — > bb^^ ^^{'~f)X 
gives the same signatures as the signal processes (both at leading and next-to-leading 
order, with X representing missing particles). As rrit should be around or greater than 175 
GeV top-antitop production will not take place at LEP II. On the contrary, at 

the NLC, tt-pairs will be copiously produced: the study of the top properties (especially 
around the threshold) is, in fact, one of the main physics goals of this machine. In order to 
quantify at leading order the importance of e~^e^ ti events (we do not consider non-top 
and single top-diagrams, as, for the purposes of this study, they can be safely neglected 
451)0, have run the code used in Refs. [^, for the choice of CM energies adopted 



here, with the decays f^^'^fi and W~ In order to have a realistic estimate 

of the cross section at threshold for tt-events at a y/s = 350 GeV NLC, we have taken 
nit = 174 GeV for the top mass (see Ref. [|^, where differences y/s — 2mt ~ 2 GeV were 
also considered). As we are interested in ti bbW^W~ events faking possible Higgs 
signals in the 66-channel, we look at the rates in the 2AMh window. We also apply the 
cuts in AMz, cos 9^+^- and pj.'*. Because of the relatively small BR of the charged vector 
bosons into muons (in fact, BR{W~ — > ^ 11%), the rates for bbW^W~ events 

with Mft5 around Mh = 90 GeV are rather small. For the non-radiative process one finds 
about 2 X 10~^ and 9 x 10~^ fb, at ^/s = 350 and 500 GeV, respectively. In presence of 
real hard photon radiation[], rates are even more suppressed compared to the signal ones: 
they are of the order 10~^ and 10~^ fb, respectively. In addition, tt-events have a quite 
large value of missing energy, i.e., -Emiss - 35 — 40 GeV (because of the neutrinos from the 
W-decajs, which escape the detectors), whereas the final states fi~^ fi~bb{'y) are in principle 
fully constrained: this should eventually allow a further reduction in the importance of 
top events. Therefore, we can neglect them in the present analysis. 

By looking at Tables I-III it is clear how the knowledge of the rates due to process 
d^) with a hard photon can be important in successfully carrying out a 6-tagging analysis, 
especially at NLC energies and if one adopts the direct reconstruction method. In fact, 
the signal and background LO rates obtained in the window AMh < 10 GeV are generally 

^^We also neglect QCD (and QED) Coulomb-like interactions between the two top quarks at threshold 
. Both because the large top mass implies that the typical spikes at ^/s « 2mt do not show dramatically 
up in the excitation curve and also because their inclusion would not change our conclusions about the 
importance of the tt-background. 

^^For the process e+e~ ti ^ bb^i'^ ii~^iX we have used a FORTRAN code produced by MadGraph and 
HELAS interfaced with routines generating the ^+i'^(7) and fi^9i_i{'-f) decays (again, only 

the tt-resonant graphs have been considered). 
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of the same order as the corrections that one gets from the M^^; NLO spectra, whereas the 
rates obtained from the M^^^ distributions are neghgible compared to those from M^j at 
LO. The NLO corrections are relatively larger at ^/s ~ 350 GeV than at ^/s ^ 200 GeV. 
The overall effect is a relative increase of the background component in the bb candidate 
sample, although only of a few percent (at all values of a/s). The cross sections in Tables I- 
III also show how an analysis that considers only radiative events e+e^ fi^fi^bb'j with 
Pj. > 1 GeV and uses exclusively the spectrum in M^^^ would be more complicated, as 
rates are at least one order of magnitude smaller than those obtained from the sum of 
LO and NLO cross sections, and because the relative excess of signal events is reduced. 
Furthermore, for the M^^^ rates separately, significances would be much smaller (see also 
Figs. 6-8). 

Figs. 5a-c shows the differential distribution in M^g for the non-radiative process 
e+e" — i> fi^fi'bB^^, for the signal (shaded) and the total background. The same sequence of 
cuts as in the Tables has been implemented here (apart from the restriction in the window 
Mf,i ^ Mh). Rates are shown for = 200 (Fig. 5a), 350 (Fig. 5b) and 500 GeV (Fig. 5c), 
with Mh = 90 GeV, plotting the histograms by bins of 2 GeV. In Fig. 6-8 we present the 
spectrum in Mf^i (a) and in M^^^ (b) in the case of the radiative process e~^e~ — > ji'^ji~bb^ 
(with pj, > 1 GeV), for the signal (shaded), and the total background. The combination 
of y/s and Mh is the same as in the previous three plots. 

The interpretation of Figs. 5a-c is quite straightforward: the signal clearly shows a 

narrow peak at Mh = 90 GeV, whereas the background (which includes the sum of the all 

contributions from crabs and deers of Fig. 1) has a broader structure. The ZZ component 

contributes with the typical Z bb peak, whereas the tail at small values of Mj,i is largely 

due to Z'y background events. The steep fall of M^^ around 120 GeV at y/s = 200 GeV is 

a kinematical effect due to the limited CM energy available at LEP II, whereas the long 

tail for Mftj > 90 GeV at the NLC is an effect due to ZZ (and also 7Z) events. The 

77-term as well as the deers and the various interferences are smaller and do not bring 

any substantial feature into the total spectrum. In general, deers are bigger than the 

77 crabs, and 6-deers dominate over /i-ones, especially at higher energies, whereas the 

total interference (of the background) is at the same level as the 6-deers. As already 

mentioned, the interference between signal and background is negligible with respect to 

^^That is, the 2^4 process, without any ISR effect. When e.m. emission from the incoming e+e"- 
hnes is included (with < 1 GeV), shapes and normahsations are generally different. However, we have 
not plotted here the corresponding curves, since in order to obtain them in the region of interest (i.e., 
AMh < 10 GeV) it is enough to rcnormalisc those in Fig. 5, according to the rates given in Tables I-III. 
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the terms discussed above and, for simplicity, it has not been included into the figures. 

The presence of a hard photon (pj- > 1 GeV) significantly modifies the lowest order 
distributions in the 66-invariant mass (especially of the Higgs process), see Figs. 6-8. By 
looking at the M^^ spectrum in the case of the signal one easily recognises the tail at small 
invariant masses due to the 67-component, whereas the 67- and //7-diagrams contribute 
to the resonance around 90 GeV. In the case of the spectrum in M^g^ it is the other way 
round. Here, the Breit-Wigner peak comes from the amplitude squared due to graphs 5 & 

6 in Fig. 2, whereas diagrams 1 & 2 and 3 & 4 appear via the tail at values of M^i^ greater 
than 90 GeV. Although not plotted in the figures, we note that the contributions from the 
//7 diagrams are comparable to those from the 67 ones at LEP II, whereas at NLC energies 
(both at 350 and 500 GeV) the 67 rates are larger. This is clearly due to the effect of the cut 
in pt on the spectrum of the photon produced in the initial state of the e'^e~ — > ii'^ii~bb'y 
process, as this latter is clearly harder at higher values of ^/s. The contribution of 67- 
diagrams is generally smaller with respect to the previous ones by at least one order of 
magnitude, whereas the interference between the three sets is completely negligible. The 
effect of the e.m. radiation on the background is not so promptly identifiable. However, 
for the distribution in M^j one easily recognises the ZZ contribution via the peak at Mz, 
which has contributions from 7Z too. The Z7 and 77 crabs have a substantially flat 
spectrum, the same for the 6-deers, whereas the /x-deers give a 0{1%) contribution to 
the Z-peak. In the case of the Mf,^^ distribution, the Z-pcak has a long tail for Mi,^ ~ 90 
GeV (such that it practically disappears), due to bby combinations in which the photon 
does not come from a Z bb'j decay. The 77 crab and contributions are small 
and flat. In general, at NLO, crabs are more than one order of magnitude bigger than 
the deers, with the ZZ and 7Z largely dominant. Also for this process the (negligible) 
interferences between signal and background have not been plotted in the figures. 

A final remark is needed if one considers that it would generally be impossible to tag 
photons too close to 6-quarks. In fact, partons give rise to jets with a finite angular 
size, such that if the photon fails within this cone it will not be distinguished from the 
other parts of the jets. Thus, its energy is counted as part of that one of the hadronic 
system associated to one of the 6-quarks and the bb invariant mass is not experimentally 
measurable for e~^e~ fi~^fi~bby events. Clearly, this frequently happens in the case of 
diagrams 5-6 in Fig. 2, 5-6 and 12-13 in Fig. 3 and 4. Occasionally, this also occurs for 
the graphs in which the photon is emitted by Icpton lines as well as by 6-lines in diagrams 

7 & 14 of Fig. 4, when it is collincar with one of the final 6's. Whereas in the first case 
such effect could well be important, we expect it to be rather small in the second one. The 
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overall effect would be that radiative events of the type (0) would look like leading events 
(|l|) and that the significance of the signals would be probably improved (as the invariant 
mass of the hadronic system embodying the untagged photon would now reproduce the 
Higgs mass, in the case of the mentioned diagrams of Fig. 2). However, in order to correctly 
predict such an effect one would need to know the terms due to the virtual corrections of 
the FSR along the 6-lines, which are not computed here. An alternative strategy, that 
we exploit, is to impose an additional cut that forces the photon to be well separated 
from both 6-quarks, for example cosOb^ < 0.95 (corresponding to a cone with an angular 
size of ~ 18 degrees). After applying this additional requirement, we get that in case of 
the signal the NLO rates are reduced by ~ 19% at LEP I, and of ~ 1% and ~ 4% at 
a NLC with ^/s = 350 and 500 GeV, respectively. The sum of the background suffers 
a reduction of ^ 7(4) [1]% at ^/s = 200(350) [500] GeV. Therefore, this approach should 
not drastically modify the significance of hard photon events, al least where these are 
quantitatively important: that is, at the NLC. 

Finally, by integrating the spectra in invariant mass at LO and NLO one obtains that 
the lowest order rates are generally increased by approximately 12(24) [27]% at ^/s = 
200(350) [500] GeV for the signal, whereas for the background from 18% (at LEP II) 
to 42(60)% (at the NLC, with ^/s = 350(500) GeV). However, as we expect that the 
corrections to the invariant mass spectra due to real photon emission behave quite similarly 
to those in Figs. 6-8 also in the case Mh ^ Mz (in the range, let us say, 60 GeV ~ Mh ~ 
120 GeV), then, when the peaks of the H and of the Z-particle are well separated (in 
the missing mass and/or in the 66-invariant mass distribution), the inclusion of radiative 
events in the sample should not change the results that one obtains by an analysis at lowest 
order (no radiation). Therefore when and Mz are not degenerate, and the difference 
SMhz = Mh — Mz is larger than approximately four times the width of the Z-boson (see, 
e.g., Ref. [^), e.m. radiative corrections can be safely neglected. On the contrary, for 
SMfjz ~ ^Tz, a full simulation is in principle needed in order to correctly establish the 
excess of bb decays in the fi~bb{'y) sample, especially at NLC energies. 

4. Summary and conclusions 

We have studied next-to-leading order electromagnetic contributions via hard photon 
radiation to the process e~^e^ fi^fi^bb, whose signature represents the 'golden channel' 
to detect and study the intermediate mass Higgs boson of the Standard Model at e+e~ 
colliders of the present (LEP II) and future (NLC) generation. In fact, such events involve 
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the production of the Higgs particle via the bremsstrahlung process e+e~ — > ZH followed 
by the decays Z — > ji^iJr and H — > hh. We restricted our attention to the ZH process 
only, thus neglecting Higgs production via the WW- and ZZ-fusion mechanisms, because 
of the particular importance of the first one in investigating the properties of the if-scalar. 
We also focused on the case of mass degeneracy Mh ~ Mz- In such conditions, a Higgs 
signal can be searched for by using 6-tagging techniques, to establish in ii'^ii~bb samples 
a content of 6-quarks much larger than the one arising from background Z ^ bb decays 
only. 

Events of the type e'^e" — > ^'^^"bb^ should be included in the phenomenological 
analysis, because at LEP H and NLC energies the initial state itself produces many hard 
photons (contrary to LEP I, where their emission is suppressed by the Z-width). Since 
both muons and 6-quarks are charged particles, photons can also be radiated in the final 
state. As it was not a priori clear how this radiation modified the integrated and differential 
rates of signal and background events and their interplay, it was important to investigate 
its effects. 

Two search strategies have been considered: a missing mass analysis of the system 
recoiling against the Z-boson identified via the muon pair, and the direct reconstruction 
of the Higgs peak via the 66-pair. Whereas for the first procedure the inclusion of hard 
photons does not modify the Miccoii spectra that one gets from the 2^4 process corrected 
for soft and coUinear ISR, we found that, in the second case, such radiative events give 
contributions comparable to the differences expected at leading order between the H ^ bb 
and the Z ^ bb rates, for ^/s ~ 350 GeV, when a window around Mh — 90 GeV is selected. 
In particular, NLO contributions are generally larger in the case of the signal, although 
the overall effect is a slight reduction of the expected relative excess of 66-pairs from the 
Higgs decay. At LEP II the influence of higher order corrections is smaller and, in first 
instance, negligible (in this case NLO rates are even the same for signal and background). 

Therefore, in an analysis that uses the direct reconstruction method to disentangle 
Higgs signals at the NLC and in the region Mh ~ Mz, the next-to-leading order rates 
should be taken into account, in order to predict correctly the amount of H and Z de- 
cays into 66-pairs. In the case of a missing mass analysis, higher order electromagnetic 
contributions can be always safely neglected. 

For values of Mh different from Mz, when the peaks of the two particles are well 
separated, the inclusion of the e . m . hard radiation should not change the conclusions one 
gets at lowest order. The only effect is an overall correction to the normalisations of the 
differential distributions, which should not drastically modify the significance factors of 
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the signal-to-background analysis. 

Before concluding, we remind the reader that our analysis did not make use of a full 
calculation up to the 0{al^) order^^ Only contributions due to the tree-level 0{a'^^) 
process e"'"e~ — > fi'^fi'bb and its ~ corrections due to the complete universal ISR (both 
real and virtual) and to real photon emission from the final state (with pj. > 1 GeV) have 
been included here. Therefore, a systematic error related to the ignorance of the virtual 
FSR effects and of the non factorisable corrections comes with our results. However, we 
notice that some of these corrections should go in the same direction as our conclusions, 
that is, of the importance of 0{aem) effects on the invariant mass spectra. For example, 
according to the Kinoshita-Lee-Nauenberg theorem |^^, one expects that logarithmic 
contributions due to 'coUinear' emission of photons from the final state must not appear 
in the expression of the cross sections, these being canceled by the (negative) contributions 
due to virtual final state photons. In the present paper, the former are included^ whereas 
the latter are not: therefore, the size of what we called NLO terms over-estimates that of 
the real ~ aem FSR corrections in the total cross section. Conversely, since such negative 
terms have the same kinematics as the 2—^4 process, they tend to reduce the LO piece. 
In particular, the Higgs peaks in the M^^ distributions should be less pronounced in the 
end. Thus, neglecting these virtual terms corresponds to somewhat under-estimating the 
effect of hard photon emission on the distributions that are of interest for Higgs searches. 
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Table Captions 

Tab. I Cross sections in femtobarns in the window |Mj,^(^-) — Mj^l < 10 GeV (being = 
90 GeV), after the cuts: pip^ > 1 GeV, |M^+^--M^| < 10 GeV and | cos^^+^-| < 0.8, 
for signal (S) and background (B), together with the total significance S (see page 
10-11 in the text for their definitions), for the integrated luminosity C = 0.5 fb~^. 
The label LO refers to the rates of the leading process (|l|) corrected for the ISR 
due to soft, collinear and virtual photons, whereas the label NLO identifies those 
from process @ with hard photons {p^ > 1 GeV). The second line for NLO events 
represents rates derived from the distributions in M^g^, whereas the others refer to 
cross sections and significances obtained from the M^^ spectra. The CM energy is 
y/s = 200 GeV. 

Tab. II Same as Tab. I, with = 350 GeV and £ = 20 fb'^ 
Tab. Ill Same as Tab. I, with = 500 GeV and C = 20 fh-\ 

Figure Captions 

Fig. 1 Feynman diagrams at tree-level for the process e^e~ fi'^^~bb (with the la- 
belling: (1), e~ (2), yU+ (3), fi~ (4), b (5) and b (6)): a) crabs; b) /i-deers; c) 
6-deers and d) the signal ZH. Internal wavy lines represent a 7 or a Z in the case of 
the graphs a, b and c, whereas for the signal they represent a Z only. Permutations 
of the virtual boson lines along the fermion ones are not shown. 

Fig. 2 Feynman diagrams at tree-level for the process e+e^ fi^fi^bb'-f (with the la- 
belling: e+ (1), e~ (2), yU+ (3), fi~ (4), b (5), b (6) and 7 (7)) in the case of the signal 
ZH. Internal wavy lines represent a Z. 

Fig. 3 Feynman diagrams at tree-level for the process e~^e~ fi'^fi~bb'-f (with the la- 
belling: e"*" (1), e" (2), fi~^ (3), fi~ (4), b (5), b (6) and 7 (7)) in the case of the crabs. 
Internal wavy lines represent a 7 or a Z. 

Fig. 4 Feynman diagrams at tree-level for the process e+e^ fi^^^bb'-f (with the la- 
belling: e"*" (1), e" (2), /i+ (3), /i~ (4), b (5), b (6) and 7 (7)) in the case of the deers. 
Labels (for i = 3,4,5 and 6) identify the 6[/i]-deer contributions. Internal wavy 
lines represent a 7 or a Z. 
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Fig. 5 Distribution in the invariant mass of the 66-pair (M^^) in non-radiative e^e — > 
jjf^jjrhh events, for signal (shaded) and background (see page 10 in the text for their 
definitions), for Mh = 90 GeV, ^ = 200 GeV (a), ^ = 350 GeV (b) and ^ = 500 
GeV (c), after the sequence of cuts: > 1 GeV, |M^+^- - Mz\ < 10 GeV and 
I cos9n+n- 1 < 0.8. 

Fig. 6 Distribution in the invariant mass of the hb- {M^ji, a) and of the 667-SYstcm (Mj^g^, 
b) in e+e^ fi~^fi~bb'~f events, for signal (shaded) and background (sec page 10 in 
the text for their definitions), for = 90 GeV and y/s = 200 GeV, after the 
sequence of cuts: pf^ > 1 GeV, |M^+^- - Mz\ < 10 GeV and | cos6'^+^-| < 0.8. 

Fig. 7 Same as Fig. 6, with ^/s = 350 GeV. 

Fig. 8 Same as Fig. 6, with = 500 GeV. 
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Fig. 2 
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Diagrams by MadGraph 
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Fig. 4 
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